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OS : Abstract. We study the narrow-line region (NLR) of the Seyfert-2 galaxy NGC 1386 by means of long-slit spectroscopy 

obtained with FORS 1 at the VLT. We use the galaxy itself for subtracting the stellar template, applying reddening corrections to 
fit the stellar template to the spectra of the NLR. The continuum gets steadily redder towards the nucleus. The spatial distribution 
of the reddening derived from the Balmer decrement differs from the continuum reddening, indicating dust within the NLR 
' with a varying column density along the line of sight. Using spatially resolved spectral diagnostics, we find a transition between 

, central line ratios falling into the AGN regime and outer ones in the H ii-region regime, occuring at a radius of ~6" (310 pc) in 

all three diagnostic diagrams. Applying CLOUDY photoionisation models, we show that the observed distinction between H ii- 
like and AGN-like ratios in NGC 1386 represents a true difference in ionisation source and cannot be explained by variations 
lO ■ of physical parameters such as ionisation parameter, electron density or metallicity. We interpret it as a real border between the 

' NLR, i.e. the central AGN-photoionised region and surrounding H ii regions. We derive surface brightness, electron density, and 

ionisation parameter as a function of distance from the nucleus. Both the electron density and the ionisation parameter decrease 
^ with radius. We discuss the consequences of these observations for the interpretation of the empirical NLR size-luminosity 

relation. In the outer part of the NLR, we find evidence for shocks, resulting in a secondary peak of the electron-density and 
ionisation-parameter distribution north of the nucleus. We compare the NLR velocity curve with the stellar one and discuss the 
differences. 
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1- Introduction covered from HST narrow-band images of seven radio-quiet 

^ PG quasars (Bennert et al., 2002). However, the slope is 

1.1. The narrow-line region in active galaxies discussed controversially: While Bennert et al. (2002) find 

The narrow-line region (NLR) in active galactic nuclei (AGNs) ^[omy suggesting a self-regulating mechanism that 

is thought to be ionised by the luminous central engine, most determines the size to scale with the ionisation parameter, 

likely an accreting supermassive black hole (BH). The narrow ^chmitt et al. (2003b) report a relation of oc L^^^^^ for their 

emission-lines are used in essentially all optical surveys to find sample of 60 Seyfert galaxies imaged with HST Such a slope 

and classify AGNs. Emission-line diagnostic provides a pow- expected for gas ionised by a central source (Stromgren) in 

erful tool to study the size and properties of the NLR. However, ^^^^ °^ constant density, 
until today, even basic issues remain open: For example, what 
determines the NLR size and structure? What are the physical 

conditions within the NLR such as electron density and ioni- S^^^^ying the [O m] emission alone as a tracer of the NLR 

.- . J J »u tt, J- t f »u size and structure has the major drawback that this emission 
sation parameter and do they vary with distance from the nu- 

^jg^^^ can be contaminated by contributions from starbursts, shock- 



Recently, a relation between NLR size and 
[Om]/15007A (hereafter [Oiii]) luminosity has been dis- 



ionised gas or tidal tails, resulting in an apparent increase of 
the NLR. In addition, different sensitivities can lead to differ- 
ent size measurements: When comparing groundbased [Oiii] 
Send offprint requests to: Nicola Bennert images of Seyfert galaxies from Mulchaey et al. (1996a) with 

* Based on observations made with ESQ Telescopes at the Cerro the HST snapshot survey of Schmitt et al. (2003a), the latter 
Paranal Observatory under programme ID 72.B-0144 reveal, on average, six times smaller NLR sizes, probably due 
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to the 15 to 20 times lower sensitivity. These considerations 
question the definition "NLR size" from [O m] imaging alone. 

Long-slit spectroscopy is a valuable alternative approach as 
it can directly probe the size in terms of AGN photoionisation 
and discriminate the contribution of stellar or shock ionisation. 
In addition, several physical properties of the NLR such as the 
electron density and the ionisation parameter can be measured 
as a function of distance from the nucleus. Thus, it allows e.g. 
to probe the constant density needed to explain the 0.33 slope 
of the NLR size-luminosity relation. These parameters are also 
of general interest as they are important input parameters of 
photoionisation models of the NLR. 

A well suited object for such an approach is the nearby (vhei 
= 868 ± 5 kms-\ NED^) Seyfert-2 galaxy NGC 1386, allow- 
ing for a detailed study of the extended NLR down to small 
spatial scales (1" ^ 52 pc). 

We here describe long-slit spectroscopic observations ob- 
tained with FORSl' at the VLT^ and the methods used to de- 
termine the "real" NLR size as well as physical conditions 
within the NLR. These methods were applied to a larger sam- 
ple of Seyfert galaxies, revealing similar results as those for 
NGC 1386 which wiU be discussed in a subsequent paper. 

1.2. NGC1386 

NGC 1386 is one of the nearest known Seyfert galaxies (D ~ 
1 1 Mpc for V3K = 774 km s"! and Hq = 71 km s'^ Mpc"'; 1" ~ 
52 pc). 

Due to its relatively high incUnation (i ^ 71°), the morpho- 
logical galaxy type has been a matter of debate: While several 
authors favour SO, an Sa classification is assumed by Weaver 
et al. (1991). Malkan et al. (1998) give a morphological type 
of Sb/c based on HST imaging, noting that they possibly have 
missed bars on larger spatial scales. In RC3 (De Vaucouleurs et 
al., 1991), NGC 1386 is indeed classified as barred SBO galaxy. 
Also NED gives SB(s)0-i- as morphological classification. 

NGC 1386 has been investigated by various authors. We 
here summarise the most important results with respect to our 
study. 

Weaver et al. (1991) study the extended NLR by means of 
groundbased imaging as well as long-slit spectroscopy with 
multiple slit positions. They show a 2D coverage of the inner 
kpc of this galaxy and discuss the overall morphology, veloc- 
ity field, electron-density distribution, and ionisation structure. 
They use the galaxy itself as stellar template to correct for un- 
derlying absorption lines. Using one diagnostic line-ratio dia- 
gram ([Oiii]/H;6 versus [Nii] /16583 A/Ho'), they determine the 
NLR (i.e. the AGN powered region) to extend ~6" to the north 
and south of the nucleus while emission lines further out can be 
attributed to Hii regions. They find an electron density which 
is decreasing with radius. The velocity field is interpreted in 
terms of a combination of a normally rotating component and 
a component undergoing high-velocity infall or outflow. 



' NASA/IPAC Extragalactic Database 

^ FOcal Reducer/low dispersion Spectrograph 

^ Very Large Telescope, Cerro Parana!, Chile (ESO) 



Mauder et al. (1992) resolve the NLR down to a linear scale 
of C.'S (^ 15 pc) using speckle interferometry and detect indi- 
vidual NLR clouds on these scales. They argue that the ionising 
radiation must be absorbed on scales <15 pc in a clumpy struc- 
ture of the NLR. 

Tsvetanov & Petrosian (1995) give positions and fluxes of 
44 H n regions distributed in a ring seen in their continuum- 
subtracted HQr-i-[Nn]/l/l6548,6583A (hereafter [Nn]) image. 
The distances of the Hii regions are ~6-12" from the nucleus 
at a position angle (p.a.) of the major axis of the ring of 25°, 
comparable to the p.a. of the major axis of the galaxy [25°; RC3 
(De Vaucouleurs et al., 1991)]. 

Ferruit et al. (2000) present HST images of the [O iii] and 
Ha+[Nn] emission. The images show the presence of very 
strong dust features, especially on the north-west side of the 
galaxy and in the nuclear regions. They find gradients in the 
[O in]/([N n] -I- Ha) ratio which at least partly cannot be due to 
dust lanes but may correspond to a real transverse change in the 
excitation conditions of the ionised gas. Ferruit et al. (2000) re- 
port a faint inclined ring of Hn-region emission extending out 
to ~12" from the nucleus with a p.a. of 25°, in agreement with 
the results from Tsvetanov & Petrosian (1995). The [Om] im- 
age of Ferruit et al. (2000) is included in the study of Schmitt 
et al. (2003a) and presented here in Fig. 1 with the slit position 
of our observation overlayed (p.a.[oiii] = 5°). It traces the di- 
rection of the maximum extent in [O m] which is displaced by 
20° from the host galaxy major axis (p.a.gaiaxy - 25°) and by 
~15° from the extended radio emission [p.a.radio = 170°, Nagar 
etal. (1999)]. 

Rossa et al. (2000) and Schulz & Henkel (2003) study the 
kinematics of the inner emission-line region in detail. Both 
find evidence for a regular velocity field showing signs of ro- 
tation. The line profiles show pronounced substructure sugges- 
tive of non-circular motions like locally expanding gas sys- 
tems. Schulz & Henkel (2003) suggest a near-edge-on warped 
rotating spiral disk as traced by Ha with a central velocity gra- 
dient corresponding to a Keplerian mass estimate of ~ 5 • 10^ 
Mq inside 0.8 kpc. Like Mauder et al. (1992), Rossa et al. 
(2000) detect several individual NLR components. 

NGC 1386 is also included in the sample of 18 Seyfert-2 
galaxies studied with long-slit spectroscopy by Fraquelli et al. 
(2003). They present reddening and surface-brightness distri- 
butions as well as a decreasing electron density observed along 
a p.a. of 169°. However, their data have a significant lower 
signal-to-noise ratio (S/N) than our data and, moreover, they 
did not take into account the underlying absorption owing to 
the contribution of the stellar population which we show to be 
important. 

2. Observations 

The high S/N long-sUt spectra of NGC 1386 were obtained us- 
ing FORSl attached to the Cassegrain focus of UTl at the VLT 
on the 25th of February 2004. Observations were made in the 
spectral range 3050-8300 A flirough flie nucleus of NGC 1386 
with exposure times of 1800 s (to study the emission Unes far 
from the nucleus) and 30 s (to gain the central emission fines 
which were saturated in the 1800 s exposure) with a typical 
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seeing of under the use of the atmospheric dispersion cor- 
rector. The detector used is a 2048 x 2048 pixels Tektronix 
CCD with 24yum wide square pixels. The spatial resolution el- 
ement is 0'.'2pix ^ The slit width corresponds to on the 
sky projecting to a spectral resolution of ~8 A (~450 kms~^) 
as is confirmed by the full-width-at-half-maximum (FWHM) 
of wavelength caUbration lines and the [O i] A5511 A night-sky 
Une. The long slit was orientated along the p. a. of the maximum 
extent of the high excitation gas observed in the [O ni] image 
(p.a.[oni] = +5°; Fig. 1) and corresponds to 6f8 in the sky. 



NdiC 1386 



1 

1 




4 2 -2 -4 

RA (arcsec) 

Fig. 1. [Om] continuum-subtracted image of NGC 1386 taken with 
HST, kindly provided by Henrique Schmitt. The contours start at the 

3iT level above the background (2.1 ■ 10"'"* ergcm"^ s"' arcsec"^) and 
increase in powers of 2 times 3cr (3cr x 2"). The exposure time is 800 s. 
The pixel size is 0'.'0455. The [Om] emission was centered on the 
optical nucleus determined from the continuum image by Schmitt et 
al. (2003a). The velocity curve suggests a hidden nucleus at ~ 1" north 
of the optical nucleus. The slit position of our observations at p.a. = 
5°, tracing the major axis of the [Om] extension, is indicated by the 
dashed line. 



3. Reduction and Analysis 

3.1. Data reduction 

Standard reduction including bias subtraction, flat-field cor- 
rection, and cosmic-ray rejection was performed using the 
ESO MIDAS^ software (version Nov. 99). Night-sky spectra 
at r-3' distance on both sides of any notable galaxy emis- 
sion were interpolated in the region of the galactic spec- 
trum and subtracted. The spectra were rebinned to a scale 
of 2.65 A pix ' during wavelength calibration. The curve of 
Tiig (1997) was used to correct for atmospheric extinction. 
The spectra were flux calibrated using the standard star CD- 

^ Munich Image Data Analysis System, trade mark of ESO 



32°9927. Foreground Milky Way reddening was corrected us- 
ing values from Schlegel et al. (1998) as listed in NED {E(b-v),g 
- 0.01 2 mag) and the extinction law from Savage et al. (1979). 
Forbidden-line wavelengths were taken from Bowen (1960). 
A heliocentric correction of -28.8 km s"^ was added to the ob- 
served velocities. 

In the spatial direction perpendicular to the dispersion axis, 
five pixel rows were averaged from the frames cleaned in this 
way according to the seeing to enhance the S/N without loosing 
any spatial information. Thus, each resulting extracted spec- 
trum of the spatially resolved emission corresponds to 1" and 
0'.'7 along and perpendicular to the slit direction, respectively. 
We choose the spectrum with the maximum intensity of the 
continuum as "photometrical center" ("zero" on the spatial 
scale). It coincides with the highest emission-line fluxes in Ha 
and [O III]. In the following, we also refer to it as "central spec- 
trum". 

3.2. Subtracting tiie stellar population 

Starlight constitutes a substantial fraction of light gathered in 
the optical spectra of AGNs, particularly in low-luminosity ob- 
jects such as LINERs and Seyfert galaxies. The emission-fine 
fluxes are afl'ected by underlying absorption lines which can be 
large especially in the case of H/3. As a consequence, starlight 
can have a great influence on the deternunation of the redden- 
ing and the AGN diagnostic line ratio H8/[0 m]. Thus, remov- 
ing the contribution of the stellar population is one of the first 
and most critical steps in the analysis of AGN emission-line 
spectra. The most widely used technique for starUght subtrac- 
tion makes use of an appropriately chosen spectrum of a normal 
galaxy. However, the stellar populations in the inner regions of 
AGNs present a variety of characteristics (Cid Femandes et al., 
1998), giving rise to the need of a library with a large diversity 
of stellar populations. Alternatively, the shape and strength of 
the starlight component can be estimated by means of stellar 
population synthesis techniques, following the work pioneered 
by Bica (1988). The classification refers to the library of syn- 
thetic spectra of Bica (1988), i.e. templates composed of dif- 
ferent percentages of star cluster spectra of several ages and 
metalficities. An S3 template, for example, contains only the 
contribution of stars older than 10 Gyr. In S4, 95% of the light 
at 5870 A comes from old stars and 5% from young and inter- 
mediate age ones. In S5, the contribution is 85% from old stars 
and 15% from young and intermediate age ones, etc. Thus, a 
low number "Sx" indicates a rather old stellar population, a 
high number a young one. This method, however, involves sev- 
eral assumptions of stellar modefing and does not necessarily 
yield an unambiguous solution. 

Here, we appUed the straightforward approach of using the 
galaxy itself as starlight template which should give the best 
representation of the stellar population. The high-sensitivity 
spectrum allows us to average several pixel rows with high S/N 
in the outer part of the galaxy spectrum. One assumption is that 
the stellar population does not vary much in the inner few ar- 
seconds, i.e. the bulge, as is confirmed by the close match of the 
stellar template and the absorption lines in aU studied spectra. 
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Furthermore, we assume that the contribution of the AGN fea- 
tureless continuum is negligible to not overestimate the under- 
lying absorption lines by scaling the template to the continuum 
in each spectrum. This seems to be justified for type-2 AGNs: 
Cid Femandes et al. (1998) find little or no dilution of the stel- 
lar lines by an underlying featureless continuum for most of the 
20 Seyfert-2s in their sample, including NGC 1386. Weaver et 
al. (1991) apply a similar method of subtracting the stellar con- 
tribution in NGC 1386. 

The stellar template was gained at a distance of 
~18" (940 pc) north from the nucleus, averaged over 5" and 
median filtered along the spectral axis over three pixels to in- 
crease the S/N. The extracted region is both far enough from 
the central region to exclude any NLR emission, as well as near 
enough to still guarantee a comparable stellar population with 
high S/N. This assumption is strenghtened by the stellar popu- 
lation analysis carried out for NGC 1386 by Cid Femandes et 
al. (1998): The equivalent widths at the nucleus are typical of 
an S3 template out to 10" where the stellar population resem- 
bles S4-S5 due to the presence of spiral arms. This coincides 
with the inclined ring of emission seen in the Ha'+[Nnl im- 
age by Ferruit et al. (2000). Further out, the stellar population 
is again S3, and we can thus use the template determined at 
~ 18" to fit the nuclear spectra out to 10". The resulting tem- 
plate spectrum is clearly dominated by absorption fines show- 
ing strong Can H+K lines as well as the Balmer series (Fig. 2, 
middle spectrum). 

The template was scaled to each row of the NLR spectrum 
by normalization in the red (~5400-5700A), justified by the 
fact that the slope at A > 5400 A does not change significantly 
for difi'erent stellar populations (Bica et al., 1986). (Note that 
we also chose this range as it does not cover any strong NLR 
emission lines.) To allow for a possible reddening difference 
of the template and each NLR spectrum due to difi^erent dust 
amounts in different galactic regions, we applied a redden- 
ing correction to the template by fitting the continuum slope 
of the template to each NLR spectrum [MIDAS command "ex- 
tinct/long" with extinction-law from Savage et al. (1979)|. This 
procedure was often necessary to avoid a template continuum 
exceeding that of the NLR continuum at some wavelength after 
normalization. 

In Fig. 3, we show the template with and without redden- 
ing correction with respect to the nucleus. The scaled (de-) 
reddened template was subtracted to gain the pure emission- 
line spectra. Finally, any remaining AGN continuum contribu- 
tion was fit interactively by a continuum slope and subtracted 
(Fig. 2, lower spectrum). 

3.3. Emission-line fluxes and reddening 

The fluxes of the pure emission-line spectra were measured as 
a function of distance from the nucleus by integrating along 
a Gaussian fit to the line profile. The fit routine "fit/spec" 
(Rousset, 1992) was used for this purpose. In the three central 
spectra (center ±1"), the strongest emission lines were satu- 
rated in the 1800 s exposure and aU fluxes were derived from 
the 30 s nuclear exposure. 
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Fig. 2. Spectrum at a distance of 4" north of the nucleus. The upper 
spectrum shows strong Can H+K absorption hues and Balmer fea- 
tures. A template spectrum of the stellar contribution derived from the 
galaxy itself at 18" north of the nucleus was subtracted (middle spec- 
trum). The difference spectrum is the lower one (with the stronger 
emission lines truncated). Both upper spectra are shifted by an arbi- 
trary amount for comparison. 




4000 4500 5000 5500 6000 



X [A] 

Fig. 3. Central spectrum of NGC 1386 with stellar template overlay ed 
as dotted line. In the upper spectra, which were shifted by an arbitrary 
amount for comparison, the reddening of the nucleus with respect to 
the template can be clearly seen. In the lower spectra, a reddening 
correction was applied to the template to improve the match of the 
template and the continuum of the nuclear spectrum. 

The uncertainties in deriving the fluxes were mostly caused 
by uncertainties in the placement of the continuum and were 
thus estimated as the product of the FWHM of the fine and 
the root-mean square deviation of the local continuum fluxes. 
Gaussian error propagation was used to calculate the errors of 
subsequent parameters such as line ratios, ionisation parame- 
ter, etc. The resulting erros are in the range of ~ 1-15%. Note 
that we did not take into account uncertainties from stellar ab- 
sorption correction and the quality of the Gaussian fits which 
was very good given the low spectral resolution of our data. 

The spectra were dereddened using the recombination 
value for the intensity ratio Ha/Hfi = 2.87 (a typical value for T 
= 10000 K, Osterbrock (1989), Table 4.2) and an average red- 
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Table 1. Observed and reddening-corrected line-intensity ratios rel- 
ative to Hfi for the nuclear spectrum. 



Line 




-^dered 


[Oil] i3727A 


1.81 ± 


0.02 


2.63 


[Neiii]i3869A 


0.77 + 


0.04 


1.07 


He+[Neiii]i3967A 


0.48 ± 


0.05 


0.64 


[0 III] ^4363 A 


0.19 ± 


0.02 


0.23 


HeiU4686A 


0.46 ± 


0.05 


0.48 


[Om]A5007A 


11.34 ± 


0.1 


10.73 


[Fevil] i5721 A 


0.29 ± 


0.04 


0.22 


[Fevil] ^6087 A 


0.44 ± 


0.07 


0.31 


[Oi]i6300A 


0.46 ± 


0.02 


0.3 


[Fex]A6375A 


0.07 + 


0.01 


0.05 


Ha 


4.7 ± 


0.05 


2.87 


[Nil] ^6583 A 


5.6 ± 


0.06 


3.41 


[S ii]i6716A 


1.04 ± 


0.05 


0.62 


[S ii]i6731 A 


1.29 + 


0.07 


0.77 



dening curve (Osterbrock (1989), Table 7.2). Note that in the 
following, only those spectra are used which have emission- 
Une fluxes exceeding the S/N ratio of 3. 

4. Results and Discussion 

While a spectrum of NGC 1386 at 4" north of the nucleus as 
well as a nuclear spectrum were already presented in Figs. 2 
and 3, Fig. 4 additionally shows representative spectra at 5" and 
8" north of the nucleus to demonstrate the emission-line vari- 
ations with distance from the nucleus. 

In Table 1, we present both the observed and reddening- 
corrected hne-intensity ratios relative to Ufi from the nu- 
clear spectrum (uncorrected for sUt losses). For pairs of lines 
([Ora], [Oi], and [Nii]) with a fixed Une ratio (~3:1), only the 
brighter line is listed. The reddening-corrected Hfi luminosity 
is(1.42±0.01)-10'*" ergs-i. 

From the [Oiii](/l4959 A-i-/l5007A)M4363 A emission- 
Une ratio, we determine the electron temperature T^, in the nu- 
clear spectrum to = 15650 ± 1500 K. 



4.1. Reddening distribution 

Two different measures of the reddening distribution were de- 
rived: (i) the above mentioned reddening of the continuum 
slope in the central parts of NGC 1386 with respect to the tem- 
plate derived in the outer parts of the galaxy; (ii) the reddening 
distribution obtained from the recombination value for the in- 
tensity ratio HafHfi. Both reddening distributions are shown in 
Fig. 5. 

As already noted by Weaver et al. (1991), the continuum 
gets steadily redder towards the nucleus. The reddest contin- 
uum is found 1" north of the nucleus. We observe a significant 
drop 3-4" north of the nucleus not reported by Weaver et al. 
(1991). It can be directly attributed to a "blue" continuum seen 
3" north of the nucleus in the [F547M/F791W] color map by 
Ferruit et al. (2000). The blue continuum is associated with the 
northern tip of the central emission-hne structure in the [O iii] 
image (Fig. 1). 



5" N (NLR) 



jt^vJ(V^>^^V^Af^Mv'|.'Y.vy^^~Hv;|^ 



4000 4500 SOOO 5500 6000 6500 

A [A] 




4000 4500 5000 5500 6000 6500 

X Ik] 

Fig. 4. Spectra of NGC 1386 at a distance of 5" north of the nucleus 
(upper panel) as well as 8" north of the nucleus (lower panel), respec- 
tively. While at 5" north of the nucleus, the lines show typical AGN 
ratios (i.e. [Oin]/H/3 > 3, [Nii]/Ha > 1), at 8" distance, the ratios are 
typical for H ii regions (i.e. [O in]/Hj8 < 3, [N ii]/Ha < 1). 



As the match between the absorption lines of the stellar 
template and those seen in the spectra is quite close, we believe 
that reddening by dust is the cause of the spatially varying con- 
tinuum slope and not an intrinsically redder stellar population 
in the central part. 

We cannot compare the absolute values of E^b-v) directly as 
the reddening determined from the continuum slope is a value 
relative to the reddening of the outer region in the galaxy, i.e. 
the template. (For comparison, the smallest reddening correc- 
tion E(B-v) needed to fit the continuum of the stellar template 
to that of the observed centrjil spectra was set arbitrarily to zero 
in Fig. 5, upper panel.) However, the reddening value obtained 
from the emission-line ratio Ha/Hj8 does not show a compa- 
rable regular distribution. Moreover, the reddening of the con- 
tinuum slope covers a range of ~ 0.3 in ^(b-v), while a three 
times larger range is covered in the reddening distribution de- 
rived from the Balmer decrement (AE^b-v) ~ !)• This can be 
due to extinction by foreground dust in e.g. the host galaxy 
which affects both the template and the central spectra and thus 
do not reflect in the relative reddening value determined from 
the difference in continuum slope of the template and the NLR 
spectra. 
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It seems that the stellar population and the NLR are suf- 
fering different dust extinctions. Possible explanations are (i) 
high density leading to an additional collisional contribution to 
the observed Ha/Hj6 ratio independent from dust. However, ex- 
tremely high densities are needed to explain the large ratios and 
there is no observational evidence for such high densities, (ii) 
Patchy dust clouds beyond the NLR may absorb the emission 
Unes along the line-of-sight of the NLR clouds. Also this ex- 
planation is rather unlikely. We favor scenario (iii): The dust is 
intrinsic to the NLR clouds and its column density varies along 
the line of sight, resulting in the observed differences in the 
reddening of the NLR and the stellar template. 

We use the reddening distribution determined from the 
HafHfi emission-line ratio to correct for the intrinsic redden- 
ing of the NLR itself as these Unes originate in the NLR and 
thus give a better estimate for the reddening within the NLR 
than the one determined from the continuum slope. 
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Fig. 5. Upper panel: Reddening distribution of the continuum of the 
central spectra relative to that of the stellar template. Lower panel: 
Reddening distribution derived from the recombination value Ka/Ufi. 
In both panels, the edge of the NLR as determined from the diagnostic 
diagrams is indicated by dotted lines. 



4.2. Spatially resolved spectral diagnostics 

Diagnostic line-ratio diagrams of the type pioneered by 
Baldwin et al. (1981) are commonly used to distinguish 
between emission-line object classes (e.g. Seyfert galaxies, 
LINERs, Starbursts, transition objects), referring to a "total" 
spectrum or the central spectrum of one object. Here, we make 
use of the spatially resolved spectra to discriminate between 
different regions in the same object, for example the NLR 
and circumnuclear or extended H ii regions. Thus, to probe the 
"real" NLR size, i.e. the central region which is photoionised 
by the AGN, and to discriminate the contribution from circum- 
nuclear starbursts, we determined spatially resolved emission- 
Une ratios. 

In Fig. 6, we present three diagnostic diagrams of 
NGC 1386. The high S/N ratio of our spectra enables us to 
measure line ratios for all three diagrams out to 10" from the 
nucleus. The symbols are chosen such that "0" refers to the 
central spectrum, the small letters mark southern regions, the 
capital ones northern regions ("a, A" - \" distance from the nu- 
cleus = 52pc, "b,B" = 2" distance = 104 pc etc.). In the second 
diagnostic diagram, the data points of the outer 7-10" are up- 
per limits, due to the faintness of the [Oi] /16300A line (here- 
after [O i]) involved. The line ratios of the central regions up 
to 6" lie in the AGN regime (0, A/a - F/f). The outer line ra- 
tios extending to 10" (G/g - J/i) all lie in the lower left corner 
usually covered by H n regions. This is a rather sharp transition 
which can be interpreted as the edge of the AGN photoionised 
region, i.e. the NLR, occuring at a radius of 6" (310 pc), while 
emission extending to 10" originates from circumnuclear Hn 
regions. These H ii regions can be attributed to the H ii regions 
seen by Tsvetanov & Petrosian (1995) in their Ha-i-[Nn] im- 
age, occuring at distances of ~6-H" from the nucleus along a 
p.a. of 5° as used in our observations. 

Our results agree with that of Weaver et al. (1991) who re- 
port an NLR extension of ~6" based on the third diagnostic 
diagram. In the first and second diagnostic diagram, their data 
only include the nuclear region due to the faintness of the in- 
volved [O i] and [S ii] lines in the outer regions. In their third di- 
agnostic diagram. Weaver et al. (1991) find some evidence for 
high [N iiJ/Hq; plus low [O iii]/HyS ratios typical for LINERs sur- 
rounding the highest ionisation gas near the nucleus. Our data 
do not show LINER-type ratios. The discrepancy between the 
two results may be due to different stellar templates used. The 
influence of the stellar template is discussed in Appendix A. 

In Fig. 7, we show the strongest reddening-corrected line 
ratios as a function of radius. 

Compared to pure imaging, using diagnostic diagrams to 
determine the NLR size has the advantage to be less sensitive 
to flux depth and to exclude contributions of circumnuclear 
starbursts. Comparing the NLR size (radius 6" =^ 310pc) de- 
termined from our diagnostic diagrams with literature values, 
we find that it is twice as large as the one found by Schmitt et 
al. (2003a) (3", Fig. 1), possibly due to the low sensitivity of 
the HST snapshot survey (short integration time of 800 s using 
the 2.4 m HST mirror compared to 1800 s with the 8 m VLT 



N. Bennett et al: The NLR in NGC 1386 



7 



mirror and a 20 times larger pixel size^). Fraquelli et al. (2003) 
detect line emission which they classify as extended NLR out 
to a distance of ~10" from the nucleus. However, our analysis 
shows that the extended emission between 6" and 10" (310 - 
520 pc) originates from circumnuclear H n regions and cannot 
be attributed to the NLR. Note that the total [Oiii] emission 
with a S/N > 3 in our spectra extends out to r ~ 12" , but only 
the central r ~ 6" can be attributed to the NLR. 

As already pointed out by Weaver et al. (1991), the sharp 
transition can be interpreted as a real physical transition be- 
tween the AGN powered NLR and the surrounding Hn re- 
gions in the disk. H n regions may be present over the entire 
emission-line region but inside a distance of 6", the AGN ion- 
isation dominates. Therefore, it is reasonable to define the ra- 
dius determined from the transition seen in the diagnostic dia- 
grams as the size of the NLR. This conclusion is further sup- 
ported by photoionisation models presented in the following 
Section. 

4.3. NLR size versus photoionisation modeling 

The determination of the NLR size depends on the definition 
of the NLR itself. We here assume that the NLR consists of gas 
photoionised by the central AGN. Thus, we are able to deter- 
mine the NLR size from the observed transition of line ratios 
in the diagnostic diagrams from the AGN regime to that typical 
for Hn regions (Fig. 6). 

We have checked whether any mechanism exists which 
drives line ratios in the diagnostic diagrams from the AGN 
regime towards and into the regime usually covered by H ii re- 
gions, despite photoionisation by an intrinsic AGN continuum 
source. This question is of interest not only for the determina- 
tion of the NLR radius, but also for galaxy classifications based 
on diagnostic diagrams in general. 

Under the assumption of a single continuum source (the 
central AGN) being responsible for the ionisation, the follow- 
ing parameters could potentially shift Une ratios from AGN- 
towards Hn-Uke in individual diagnostic diagrams: (i) extinc- 
tion; (ii) ionisation parameter in combination with (iii) metal- 
licity gradients; or (iv) electron density. We discuss each of 
them in turn. 

(i) Extinction can be excluded to have a large effect on 
the observed Une ratios in the diagnostic diagrams. The line 
ratios in all three diagnostic diagrams were chosen to mini- 
mize reddening eff'ects by using neighbouring lines. In addi- 
tion, we measured the reddening distribution directly (Fig. 5, 
lower panel) and use these results to correct for the redden- 
ing. Using the diff'erent reddening distribution determined by 
the continuum slope (Fig. 5, upper panel) does not change the 
results within the errors. 

Models (ii)-(iv) presented below are based on the photoion- 
isation code CLOUDY [e.g. Ferland et al. (1998)] and assume a 
pointlike AGN continuum source which illuminates clouds of 
constant density. The spectral-energy distribution is a typical 
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' Note, however, that the HST observations were an imaging cam- 
paign while we carried out spectroscopy, thus these values are not di- 
rectly comparable. 
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Fig. 6. Diagnostic diagrams for spatially resolved Une ratios in 
NGC 1386. The dividing lines were taken from the analytic AGN diag- 
nostics of Kewley et al. (2001). The symbols are chosen such that "0" 
refers to the central spectrum, the small letters mark southern regions, 
the capital ones northern regions ("a,A" = 1" distance from the nucleus 
= 52 pc, "b,B" = 2" distance = 104 pc etc.). In all three diagrams, only 
the central r ~ 6" (310pc; A,a - F,f) show ratios expected for AGN- 
photoionised gas (upper right comer). Futther out (until r ~ 10"; G,g 
- J,i), the Une ratios fall into the H n-region regime (lower left comer). 



mean Seyfert continuum composed of piecewise power laws 
(Komossa & Schulz, 1997) with, in particular, an energy index 
Ouv-x = - 1 .4 in the EUV and a photon index = - 1 .9 in X- 
rays. The clouds are assumed to be ionisation bounded and of 
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Fig. 7. Reddening corrected emission-line ratios used in the diag- 
nostic diagrams as a function of distance from the nucleus: [Nii]/Hq' 
and [S ii]/Ha (upper panel) as well as [O iii]/H/3 and [O il/Ha {lower 
panel). Note that the [O iii]/H/j' ratio was divided by a factor of 30 for 
comparison. The edge of the NLR as determined from the diagnostic 
diagrams is indicated by dotted lines. 



solar metaUicity, unless stated otherwise. In Fig. 8, we present 
the effects of varying ionisation parameter (solid line 1), vary- 
ing metal abundances (dotted lines 2 & 3), varying nitrogen 
(N) and sulphur (S) abundances (dash-dotted lines 4 & 5), and 
high density (dashed line 6) on the emission-line ratios in the 
diagnostic diagrams. 

(ii) The ionisation parameter U is defined as f/ = 
Q/iAncrter^) with Q = rate of H-ionising photons, rie - elec- 
tron density, and r - distance between photoionising source 
and emission-line clouds. Variations in ionisation parameter 
alone have been studied by many authors and generally lead 
into the LINER-region of diagnostic diagrams for small ion- 
isation parameters under the assumption of solar abundances 
[e.g. Stasinska (1984); Osterbrock (1989)]. This trend is shown 
in Fig. 8 (solid line 1) which includes the range of ionisation 
parameters we measure for NGC 1386, varying from log U = 
-4.0 ... - 1.5 in steps of 0.5 (marked by circles) from bottom 
to top. Sequences were calculated three times, using densities n 
= 800, 400, and 200 cm'^ and r = 50, 100, and 250 pc, respec- 
tively (printed as circles from right to left; only the data points 



for n = 800 cm ^ and r = 50pc are connected by a line). We 
come back to (ii) in combination with (iiib) below. 

(iii) The strongest hint argueing against metalUcity varia- 
tions to change the line ratios in the observed direction is the 
similarity of all three diagnostic diagrams. In all three dia- 
grams, the line ratios of each spectrum fall in the same regime, 
either AGN or H ii region. In fact, while extreme parameters 
may shift line ratios towards H n-like in individual diagrams, 
no mechanism does so for all three, as we wiU show in the fol- 
lowing. 

(iiia) Metal abundances were varied between Z - (0.05- 
3.0) X solar. Systematically decreasing abundances first lead to 
an increase in line ratios which can be partly traced back to 
oxygen being a strong coolant, and thus an increase of heating 
in case of oxygen depletion (Fig. 8). 

The dotted lines 2 & 3 show models of varying metal abun- 
dances (3, 2.5, 1.7, 1.3, 1.0, 0.9, 0.5, 0.3, 0.1, 0.05 x solar; from 
right to left) for constant ionisation parameter. The lower dot- 
ted line 2 corresponds to log U = -2.8, the upper dotted line 3 
tologf/ = -3.7. 

(iiib) If, instead, only N and S abundances are varied, line 
ratios shift horizontally in the [N ii] and [S ii] diagrams towards 
the H II regime (Fig. 8). The dash-dotted lines 4 & 5 show mod- 
els of varying N and S abundances (3, 2.5, 1.7, 1.3, 1.0, 0.9, 
0.5, 0.3, 0.1, 0.05 X solar; from right to left) for constant ioni- 
sation parameter. The upper dash-dotted line 4 corresponds to 
log U = -2.8, the lower dash-dotted line 5 to log U = -3.7. 

Thus, a combination of outwards decreasing ionisation pa- 
rameter and decreasing metal abundances would place line ra- 
tios in the Hn-regime in the observed way, despite AGN- 
intrinsic excitation. However, the same trend does not hold for 
the [O i] diagram. Just based on oxygen lines, this diagram is 
rather insensitive to variations in N and S abundances and line 
ratios remain well within the AGN regime for the whole pa- 
rameter range (Fig. 8). We conclude that combined ionisation 
parameter and metallicity effects cannot explain the observa- 
tions. 

(iv) Finally, the presence of a high density component will 
again shift [N n] and [S ii] towards the H ii regime. To demon- 
strate an extreme, we plot in Fig. 8 the case for lognn = 6 
(dashed line 6). The ionisation parameter varies from bottom 
to top in steps of 0.5 from log U = -4.0 ... - 1.5. Again, while 
[N n] and [S ii] decrease, [O i] is strongly boosted, moving line 
ratios away from H n-like excitation. This effect can be traced 
back to the different critical densities. 

We conclude that the observed distinction between H ii-like 
and AGN-like line ratios of NGC 1386 represents a true differ- 
ence in ionisation source, and thus our method to measure the 
NLR radius is valid. The second diagnostic diagram including 
the [O i] emission-line, new in our study compared to literature 
data, was essential to reach this conclusion, since our photoion- 
isation calculations show that there are means to reach the H ii- 
regime with an AGN as ionising source in the [N n] and [S n] 
diagrams. 
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Fig. 8. CLOUDY modeling results demonstrating the effects of varying 
ionisation parameter (solid line 1), varying metal abundances (dotted 
lines 2 & 3), varying N and S abundances (dash-dotted lines 4 & 5), 
and high density (dashed line 6) (see text for details). The bold dashed 
line indicates the separation between AGN regime and Hn regions. 
Four representative line ratios of NGC 1386 are also shown, falling in 
the AGN (A, c) and Hn-region regime (J, g). 

4.4. Surface-brightness distribution 

In Fig. 9, we present the surface-brightness distribution of the 
[Om] and Ha emission line as well as of the continuum (at 
5450-5700 A). The surface-brightness distributions are similar 
to each other, centrally peaked and all decreasing with distance 



from the nucleus. They reveal a secondary peak ~3" north of 
the nucleus and a slightly lower tertiary peak ~4" south of the 
nucleus. While the [Om] surface brightness is exceeding that 
of Ha in the central parts, this behaviour changes at the edge 
of the NLR as determined from our 2D diagnostic diagrams 
(indicated by the dotted lines). 

For comparison, the [O m] surface-brightness distribution 
from the HST image of Schmitt et al. (2003a) is also shown. 
It was derived by averaging three vectorplots along the ma- 
jor axis of the NLR emission along p.a. - 5°+l°. It clearly 
shows the higher spatial resolution of the HST image (0'.'0455 
pix"') compared to the 1" spatial sampling of our spectral 
data. However, as already mentioned in Section 4.2, it once 
again reveals the low sensitivity of the HST image compared 
to our spectroscopy as it detects [O iii] emission at a S/N of 3 
only out to 3-4" from the nucleus. The HST [Om] surface- 
brightness distribution reveals several subpeaks of possibly in- 
dividual NLR clouds, as can be already seen in the [O m] image 
(Fig. 1). These substructures are smoothed out in the ~20 times 
lower spatial resolution of our spectra. 
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Fig. 9. Surface-brightness distribution of NGC 1386 in [Om] (open 
diamonds), Ho- (filled diamonds), and continuum (at 5450-5700 A, 
stars). Error bars are smaller than the symbol size. The [O m] surface- 
brightness distribution from the HST image with S/N > 3 is shown as 
small open squares connected by a line. The HST image has a 20 times 
higher spatial resolution but a significantly lower sensitivity, not al- 
lowing to measure the outer parts of the NLR. The edge of the NLR as 
determined from the diagnostic diagrams is indicated by dotted lines. 



4.5. Electron-density distribution 

Applying the classical methods outlined in Osterbrock (1989), 
we derive the electron density as a function of distance from the 
nucleus using the ratio of the [S ii] /1/167 16,6731 A (hereafter 
[Sn]) pair of emission lines (Fig. 10). We used the observed 
central temperature to correct for the dependency of electron 
density on temperature*" (Te - 15650 K). Due to the faintness 
of the involved [O iii] /14363 A emission Une, we were not able 
to measure the temperature in the outer parts. 

* n,(T) = n,(obs) • V(r/ 10000) 
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The electron density is highest at the nucleus with - 
1540 cm"^ and decreases outwards down to the low-density 
limit (assumed to be 50 cm"^). A secondary peak can be seen 
at a distance of 6" to the north of the nucleus {iie ^ 670 cmr^)? 

Comparing the position of the observed northern n^. peak 
with the 2D electron-density distribution of Weaver et al. 
(1991), the closest data point is at ~6'.'5 north-east from the 
center at roughly a p.a. of 3° with Ug < 100 cm"^ (their Fig. 
lOd). The spatial distance between this data point and our 
peak is ~25 pc. Thus, the He peak observed in our spectra may 
be attributed to an individual NLR cloud or local density inho- 
mogeneities which our line-of-sight happens to pass. Mauder 
et al. (1992) also detect individual NLR clouds using speckle 
interferometry and conclude that the ionising radiation must 
be absorbed on scales <15 pc in a clumpy structure of the 
NLR. Interestingly, the peak occurs at the edge of the NLR. 
We discuss the origin of the northern peak in greater detail in 
Section 4.7. Another, but significantly smaller peak is visible at 
the southern edge («« ^ 210 cm"^). 

While a detailed fit of a power-law function is complicated 
by the intrinsic scatter of the data points, a general dependency 
on r can be estimated. The best fits of a power-law function 
with Ueir) - riefir^^ were derived with 6 between 1 and 2 for 
the electron density (Fig. 10), neglecting the data point 6" north 
of the nucleus which may have another origin (Section 4.7). 



Electron density 
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Fig. 10. Electron density obtained from the [S n] ^67 1 6 A/A673 1 A ra- 
tio as a function of distance from the nucleus. Open symbols indicate 
locations where is in the low-density limit (assumed <50cm"'). 
The electron density decreases with radius and shows a secondary 
peak 6" north of the nucleus as well as a tertiary fainter one at 
6" south. The fits ^(.(r) = n^Qr^^ and ^(.(r) = n^QV^^ are shown as 
dashed lines. The edge of the NLR as determined from the diagnostic 
diagrams is indicated by dotted lines. 



^ Note that the temperature can be a function of distance from the 
central AGN and is most probably decreasing, if the AGN is the only 
heating source. In such a case, correcting with the central temperature 
overestimates the electron density in the outer parts. Thus, a decreas- 
ing temperature can steepen the slope of the electron density. 



4.6. lonisation-parameter distribution 

To estimate the value of the ionisation parameter, the line ra- 
tio [O n]/l3727 A/[0 m] A5007 A (hereafter [O n]/[0 m]) can be 
used [e.g. Penston et al. (1990)]. Komossa & Schulz (1997) 
have shown that in case of low density (hh < 10^ cm"^), this 
line ratio yields U independent of the shape of the ionising con- 
tinuum. It gives a good approximation of the ionisation parame- 
ter as long as there are no strong density inhomogeneities with 
radius. If the NLR contains a mixture of densities, the ratio 
predicts a too high value of U. But even if the absolute value 
of U was overpredicted, the observed slope would not change. 
We used the theoretical variations of [O n]/[0 m] calculated by 
Komossa & Schulz (1997) (their Fig. 7) with U for a power- 
law continuum (qtuv-x - -1.5) and a density of 2 < logriH 
(cm~^) < 3 (Fig. 1 1). As the electron density is decreasing with 
distance from the nucleus (~1500 to ~50cm"-'. Fig. 10), the 
resulting ionisation parameter lies in between the two curves 
calculated for logn^ = 2 and logn^ = 3, i.e. the slope nar- 
rows. 

While Weaver et al. (1991) propose a decreasing ionisation 
parameter to explain the observed line ratios in NGC 1386, we 
can determine the ionisation parameter directly and indeed find 
that it is decreasing with distance from the nucleus. It varies 
between f/iog(ne)=3 = (2.83± 0.01) -lO'^ and 2.4 • lO"'' within 
the NLR with the highest value in the center. It reveals a sec- 
ondary peak at 4" to the north of the center which is at the 
north edge of a "blob" of emission clearly visible in the [O m] 
image from Schmitt et al. (2003a) (Fig. 1). Compared to the 
electron density, the secondary peak is 2" closer to the nucleus 
and significantly broader. We discuss the origin in Section 4.7. 
In Appendix B, we discuss any possible influence of the red- 
dening distribution on the slope of the ionisation parameter. 

We estimated the general dependency of the decrease of U 
on r. The best fit of a power-law function with U(r) = Uor^'' 
is derived with 6 ~ I (Fig. 11), again neglecting the "north- 
em peak". Inserting the fit rig oc to the definition of U = 
QliAjicnei'^) yields U oc r"'. This is in agreement with the fit 
to the ionisation parameter. It indicates that the excitation of the 
bulk of the NLR is due to photoionisation and that ionisation 
by e.g. shocks is negligible. 

4.7. Northern peak: evidence for shocks 

Our observations show that both the electron density and the 
ionisation parameter are decreasing with radius. They show a 
secondary peak ~ 6" and 4" north of the center, respectively. 
The peak at 4" is close to the northern tip of the [O m] emission 
(Fig. 1). In the same region, a blue continuum indicating low 
extinction is observed [Fig. 5, upper panel, see also Ferruit et 
al. (2000)]. 

Can both the increased density and the ionisation param- 
eter originate from a shock of a radio jet? Extended radio 
emission has been found by Nagar et al. (1999), displaced 
to the observed slit p.a. tracing the major [O ni] extension by 
~15° (p.a.radio = 170° versus p.a.[oiii] =5°). A scenario in which 
the outflowing material causing the radio emission compresses 
the line-emitting gas by shocks, increasing the ambient density. 



N. Bennert et al.: The NLR in NGC 1386 



11 




0.002 - log n=2 <> 
- log n^3 ♦ 



Fig. 11. lonisation parameter derived from [O ii]/[0 iii] ratio as a func- 
tion of distance from the nucleus (open symbols: Hh = 100 cm"^, filled 
ones: Hh = 1000 cm"^). The lonisation parameter is highest in the cen- 
ter and decreases with distance with a secondary peak at 4" to the 
north of the center. The fit U(r) = J/o'' ' is shown as dashed line. 
The edge of the NLR as determined from the diagnostic diagrams is 
indicated by dotted lines. 



seems to be rather common for Seyfert galaxies [e.g. Capetti et 
al. (1996); Falcke et al. (1998)]. While the shock front results 
in a locally increased density, the continuum of the shock can 
possibly ionise a larger area of the surrounding medium and 
produce a change in the observed line ratios, resulting in an U 
bump which is broader than the sharp peak of rie. Fast autoion- 
ising shocks can produce such a local continuum ionising the 
surrounding medium (Dopita & Sutherland, 1996). 

To check the shock scenario, we searched for high- 
ionisation emission Unes (e.g. [Fex] /16375 A). In the 
nuclear spectra, we find strong [Fe vn] /15721 A and 
[Fe vn] /16087 A emission, comparable to that of [Oi]. 
Due to our low spectral resolution, [Fe x] A6375 A is 
blended with the [O i] /16363 A line. Thus, to estimate the 
contribution of [Fe x] /16375 A, we fit two Gaussians to 
[Oi] /1/16300,6363 A with line fluxes fixed to the ratio of 1:3 
and a third one to fit any remaining flux of [Fe x] A6375 A. 
Indeed, we find evidence for [Fe x] contribution in the central 
2" ([Fex]/H/3 ~ 0.047 in the center, [Fe x]/H/3 ~ 0.049 at 1" to 
the north) and at 6" ([Fe x]/HyS ~ 0.049), i.e. at the same region 
where we observe an increased electron density. 

The observed minimum in the reddening of the continuum 
slope roughly 3-4" north of the nucleus (Fig. 5, upper panel) 
could be mimicked by an additional local continuum contri- 
bution. The extinction in the NLR measured from the Ha/H/3 
value reaches a minimum at 6-7" north of the nucleus (Fig. 5, 
lower panel). This can be explained by the destruction of dust 
by local shocks. 

To conclude, we find several signs for shocks. A radio jet 
interacting with the NLR gas is a plausible explanation for 
the northern peak observed in the electron density and the 
ionisation-parameter distribution. 



4.8. NLR versus stellar velocity curve 

We derived the NLR velocity curve from the average of the 
peak wavelengths of the Ha and [Nn] emission lines. In ad- 
dition, given the high S/N ratio of our spectra, we are able 
to trace the velocity field from the "peak wavelengths" of the 
stellar absorption Une Can K (before subtraction of the stel- 
lar template) throughout the whole region as this line is not 
blended with emission lines. The uncertainty in velocity val- 
ues obtained by comparing repeated measurements at the same 
point is 20kms"'. Note that gas velocities derived from the 
emission of [O n], Hj8, and [O m] reveal a curve similar to that 
of Hot, [Nn]. 

Both, the stellar and the gaseous velocity curves, have 
a similar slope (Fig. 12). The stellar velocity curve covers 
an overall smaller range of (~720-940)±20 kms"\ while 
the emission lines show velocities ranging from (~690- 
990)±20kms"'. We define the kinematical center by taking 
the symmetry center of the outer portions of the velocity 
curve. The kinematical center has an heliocentric velocity of 
~ 825+20 km s"\ i.e. slightly lower than the one given in NED 
(vhei = 868 km s"^). However, the velocity measured at the po- 
sition of the optical nucleus is with Vhei = 870±20kms~' the 
same as the literature value within the errors. 

The kinematical center of the stellar velocity curve lies 
~1" north of the optical nucleus ("0" on the spatial scale). This 
off'set is in agreement with the analysis of the Ha peak veloc- 
ity field by Weaver et al. (1991) and Schulz & Henkel (2003) 
who report an offset of the kinematical center of ~1" to the 
north-east of the photometrical center. The kinematical center 
is suggestive of a hidden nucleus residing 1" north of the op- 
tical nucleus as has also been suggested by Rossa et al. (2000) 
(see Fig. 1). 

However, the kinematical center of the NLR velocity curve 
we find at a p. a. of 5° [compared to 23° of Schulz & Henkel 
(2003)] is ~0'.'5 south of the optical nucleus. It may be indica- 
tive of a disturbed velocity field from outflows and shocks at 
our observed slit orientation. 

While the kinematics of the emission-line region in 
NGC 1386 has been investigated by several authors, an inter- 
pretation of the observed velocity field is not unambiguous. 
Though the natural explanation is that of a slightly projected 
rotation curve from an inclined emission-line disk in the cen- 
ter, an interpretation by a bipolar flow accidentally aligned with 
the kinematical major axis of the rotating disk is also possible 
(Schulz & Henkel, 2003). However, as the velocity curves of 
both NLR gas and stars are similar, a bipolar outflow of the 
NLR is rather unlikely as it had to be decoupled from the stel- 
lar gravitational field and its accordingly stratified intersteUar 
medium. 

Instead, the NLR seems to foUow the steUar rotation in- 
dicating that the NLR gas is distributed in a disk rather than 
a sphere. Galactic rotation of narrow-line gas seems not to 
be uncommon [e.g. Schulz & Henkel (2003)]. Mulchaey et al. 
(1996b) propose that the NLR gas is distributed in the galactic 
disk rather than a sphere, to explain differences of the observed 
NLR structure with that expected from the Unified Model in 
the case of a spherical distribution of the NLR gas. 
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The gaseous velocity curve does not show a different be- 
haviour in the NLR and the H n-region regime. This indicates 
that the NLR does not consist of an intrinsically different gas 
component with its own kinematics, but is ambient gas in the 
galactic disk photoionised by the AGN. 

Our results reveal some differences between the stellar and 
the NLR velocities in NGC 1386. The NLR velocities are in 
general higher than those of the stellar absorption lines. The 
difference is especially large in both the outer southern and 
northern part of the NLR. These can be due to many effects, 
e.g. rotating gas disk/ring within a kinematically hotter stel- 
lar component (the interpretation being complicated by the fact 
that due to the high inclination of the galaxy, the hue of sight 
towards NLR intersects both the stellar disk and the bulge), gas 
moving along elliptical streamlines in a barred potential (keep 
in mind the ambiguous Hubble-type classification. Sect. 1.2), 
or outflow/jet interaction. However, for a detailed interpretation 
of the stellar and gaseous velocity field, kinematical modelling 
is required which lies beyond the scope of this paper. 
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Fig. 12. Velocity curve of the NLR derived from the average value 
of the peak wavelengths of the Ha and [Nn] emission lines (filled 
symbols) as well as the stellar velocity curve measured using the Ca ii 
K absorption-line "peak" wavelength as seen in the "raw" spectrum 
(open symbols). The edge of the NLR as determined from the diag- 
nostic diagrams is indicated by dotted lines. 



4.9. NLR size-luminosity relation 

The NLR size and its relation to [Om] luminosity has been 
studied by Bennert et al. (2002) and Schmitt et al. (2003b) 
based on [Om] images. While Bennert et al. (2002) find a re- 
lation between NLR size and luminosity with a slope of 0.5, 
Schmitt et al. (2003b) report a slope of 0.33. If the NLRs were 
ambient circunmuclear gas fully photoionised by the nuclear 
continuum source, a Stromgren-type radius-luminosity rela- 
tionship R cc /,0-33^-o.66 gjjj g slopc of 0.33 is indeed 
expected in the case of a constant density. However, a de- 
creasing density as observed for NGC 1386 seems to be com- 
mon to several Seyfert galaxies [e.g. Storchi-Bergmann et al. 



(1992); Fraquelh et al. (2000, 2003)], contradicting the simple 
"Stromgren explanation". 

A slope of 0.5 found by Bennert et al. (2002) for the NLR 
size-luminosity relation of both the [O m] emission line as well 
as Ufi is expected if on average, all AGNs have the same ion- 
isation parameter, density, and ionising spectral energy distri- 
bution (Netzer, 1990; Netzer et al., 2004). It suggests a self- 
regulating mechanism that determines the size of the NLR to 
scale rather with the ionisation parameter U: U cc Q/(R^n) 
with Q oc L(U/3) which leads to R cc L(HyS)'^^, for a given 
U and He, in agreement with the empirical relation (assum- 
ing that L(Hj6) oc L([OIII])). At the outskirts of the NLR, 
He ~ 10^-3 cm ^ and U ~ 10 ^^-^^ so that efficient [Om] 
emission comes from regions with Urie ~ 1. The NLR size- 
luminosity relation indicates that all objects have a similar line 
value at the edge of their NLR, independent of the individual 
spatial behaviour of U and Thus, the observed spatially de- 
creasing ionisation parameter and density in NGC 1386 do not 
contradict the explanation of the NLR size-luminosity relation 
in terms of a constant (in all objects) ionisation parameter. 

However, the NLR size determination using [O m] images 
alone bears various uncertainties due to e.g. the sensitivity de- 
pendency of the observation and possible contributions from 
circumnuclear starbursts as has been shown for NGC 1386. 
Compared to the spatially resolved spectral diagnostics mea- 
suring the "real" NLR size, the apparent NLR size determined 
by [O m] images can be either smaller in case of low sensi- 
tivity or larger in case of contributions of circumnuclear star- 
bursts. Using the methods described here to determine the NLR 
sizes of a larger sample of objects will help to scrutinise the 
NLR size-luminosity relation. Nevertheless, individually dif- 
fering NLR sizes may not change the overall slope of the size- 
luminosity relation which extends over four orders of magni- 
tudes in luminosity and two orders of magnitudes in size. 

5. Conclusions 

Applying diagnostic diagrams to our spatially resolved spec- 
tra, we observe a transition of emission-line ratios from the 
central AGN region to the surrounding H ii regions. This tran- 
sition occurs at r ~ 6" (310 pc) to the north and south of the 
optical nucleus and is, for the first time, observed in all three 
diagnostic diagrams, i.e. including the second diagnostic dia- 
gram involving the [O i] emission line. The most probable ex- 
planation for this transition is that the stellar ionisation field 
dominates that of the AGN at a distance of >6". We can ex- 
clude other effects such as the variation of physical parame- 
ters from CLOUDY photoionisation modelling. The second di- 
agnostic diagram ([O m]/HjS versus [O i]/Ha) was essential to 
draw this conclusion as our photoionisation calculations show 
that there are means to reach the H n-region regime with an 
AGN as ionising source in the [N n] and [S n] diagrams. We 
determine the radius of the NLR in NGC 1386 to 6" (310pc), 
independent of sensitivity and excluding [O m] contamination 
from circumnuclear starbursts. Applying this method to a larger 
sample of Seyfert galaxies will help to scrutinise the NLR size- 
luminosity relation. 
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We derive physical parameters of the NLR in NGC 1386 
such as reddening, surface brightness, electron density, and 
ionisation parameter as a function of distance from the nu- 
cleus. The differences between the reddening distributions de- 
termined from the continuum slope and the Bahner decrement 
argue in favour of dust intrinsic to the NLR clouds with varying 
column density along the line of sight. Both electron density 
and ionisation parameter decrease with radius. The decreasing 
electron density argues against a simple Stromgren behaviour 
of the NLR gas which has been suggested by the NLR size- 
luminosity relation with a slope of 0.33. The decreasing ionisa- 
tion parameter in this individual object does not rule out a com- 
mon effective ionisation parameter in AGNs as suggested by a 
slope of 0.5 of the NLR size-luminosity relation. In the outer 
part of the NLR, the secondary peak in the electron density 
and ionisation-parameter distribution are interpreted as signs 
of shocks, a scenario also supported by other observations. 

The NLR and stellar velocity fields are similar and indicate 
that the NLR gas is distributed in a disk rather than a sphere. 
The differences between the two velocity fields can be due to 
several effects such as a rotating gas disk or ring, gas moving 
along elliptical streamlines in a barred potential or outflowing 
gas from a jet interaction. 

The methods presented here were applied to a larger sam- 
ple of Seyfert galaxies, revealing similar results to that of 
NGC1386. A detailed discussion can be found in Bermert 
(2005) and wiU be summarised in an upcoming paper. 
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Appendix A: Influence of stellar template 

To probe the influence of the stellar template used, we carried 
out two approaches of (i) removing the weak emission lines of 
[On], [Oiii], and [Nii] still visible in the template and (ii) not 
using any template at aU. 

If we remove the weak emission hues in the stellar tem- 
plate (method i), the resulting Une ratios are the same within 
the errors for most spectra, i.e. the weak emission is negligi- 
ble against the strong NLR emission. The exceptions are the 
line ratios of two outer spectra (H, I) which are shifted towards 
higher [N n]/Ha values and now he at the edge between H ii 
region and LINER region. This is a possible explanation for 
the discrepancy between our results and that of Weaver et al. 
(1991) who find LINER-type line ratios in the outer region. 
They modified their stellar template derived from the galaxy 
by removing the remaining weak [N ii] and Har emission and 



replacing Ha by an artificial absorption line. This is a more 
physical approach than method (i), i.e. not only removing the 
remaining weak emission, but at the same time deepening the 
HjS and Ha absorption trough. The true absorption will prob- 
ably be underestimated due to weak emission in these Unes. 
However, such a correction is difficult to estimate and the Ha 
absorption used by Weaver et al. (1991) probably still underes- 
timates the true absorption. 

We decided to not remove the weak emission: On the one 
hand, increasing the observed [O iii] and [N n] emission by clip- 
ping the weak emission in the template and at the same time 
increasing the observed Ha and H/3 emission by deepening the 
absorption lines in the template will result in approximately 
the same line ratios and thus wiU not have a great influence on 
the diagnostic diagrams. On the other hand, the weak emission 
lines are attributable to ionised regions around young stars in 
the galaxy lying along our line-of-sight, and do not originate 
from the NLR emission. Thus, they truly contaminate the NLR 
emission and need to be subtracted. 

Applying no correction for underlying absorption at all 
(method ii) yields extremely large [O ni]/Hy6 values (>25) near 
the nucleus as has already been pointed out by Weaver et al. 
(1991). In the outer parts, the absorption trough is clearly visi- 
ble with a small H/3 emission peak in the center. In addition, the 
spectra show galaxy absorption features which accompany the 
absorption in I^. These findings show the need for correcting 
the stellar absorption Unes using a suited stellar template. 

Appendix B: Influence of reddening 

The only physical parameter in our study which is strongly sen- 
sitive to reddening is the ionisation parameter, derived using the 
line ratio [On]/[Oiiil. We here briefly discuss the influence of 
the different reddening measures on the ionisation parameter. 

Although we corrected for the reddening within the NLR 
using the recombination value Ha/Hy6, there is a clear differ- 
ence between the reddening of the continuum slope and that 
determined from the Balmer decrement: While the reddening 
distribution determined by the ratio Ha/HJS shows a rather ran- 
dom distribution, the continuum slope reddening steadily in- 
creases towards the nucleus with the exception of the blue con- 
tinuum at 3-4" to the north of the nucleus (Fig. 5). Can these 
differences explain the observed slope of the ionisation param- 
eter, i.e. the general decrease with radius and the secondary 
peak at that position at which we see a blue continuum? Is it 
possible that we do not correct for the reddening intrinsic to 
the NLR gas by using the reddening value determined by re- 
combination lines but would the continuum slope be a better 
tracer for the reddening distribution? To probe the difference, 
we used the reddening determined from the continuum slope to 
correct the observed [On]/[Oni] emission-line ratio. The re- 
sulting ionisation-parameter distribution changes only slightly 
showing an even steeper increase towards the center and a more 
pronounced secondary peak (Fig. B.l). 

Thus, the observed slope in Fig. 1 1 is not an artefact due to 
a "wrong" reddening correction from the recombination value: 
Using the reddening determined from the continuum slope, the 
ionisation-parameter distribution does not flatten but, in con- 
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trary, steepen. The general conclusion of a decreasing loni- 
sation parameter with radius which shows a secondary peak 
3^" north of the nucleus remains unchanged. However, we 
consider the recombination value Ha/Hfi as a better indicator 
of the reddening distribution within the NLR as it uses emis- 
sion lines originating in the NLR itself. We suggest that the 
differences in the reddening measurements originate from dust 
within the NLR with varying column density along the Une of 
sight. Although in general using the HafH/S value for redden- 
ing correction bears the problem of the uncertain amount of the 
underlying stellar absorption, we believe that our correction us- 
ing the galaxy itself as stellar template gives a good approach 
of the absorption-line free ratio. 
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Fig.B.l. lonisation parameter derived from [On]/[Oin] ratio as a 
function of distance from the nucleus for hh = 100 cm"^. We com- 
pare the ionisation parameter derived using the reddening determined 
from the continuum slope for correction (open symbols) with the one 
from Fig. 1 1 using the Balmer decrement to correct for the reddening 
(filled symbols). The edge of the NLR as determined from the diag- 
nostic diagrams is indicated by dotted lines. 
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